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Abstract 

The properties of neutron stars are investigated by including 8 
meson field in the Lagrangian density of modified quark-meson 
coupling model. The £~ population with 8 meson is larger than 
that without 8 meson at the beginning, but it becomes smaller 
than that without 8 meson as the appearance of S - . The 8 meson 
has opposite effects on hadronic matter with or without hyperons: 
it softens the EOSes of hadronic matter with hyperons, while it 
stiffens the EOSes of pure nucleonic matter. Furthermore, the 
leptons and the hyperons have the similar influence on 8 meson 
effects. The 8 meson increases the maximum masses of neutron 
stars. The influence of (a*,4>) on the 8 meson effects are also 
investigated. 

1. INTRODUCTION 

8 meson is an isovector scalar meson, its contribution is expected to be neglectable in nuclei with small isospin 
asymmetry and in nuclear matter at saturation density. However, for strongly isospin- asymmetric matter at high 
densities in neutron stars the contribution of the 8 field should be considered |l|. In the last decade, the effects of 
coupling to the 8 meson like field on nuclear structure properties of the drip-line nuclei, on the dynamic situations of 
■ heavy ion collisions and on asymmetric nuclear matter are investigated @, 0, 0, Q . Recently, the density dependent 
coupling constants are introduced additionally to reexamine the 8 meson effects on properties of finite nuclei and 
asymmetric nuclear matter in the Quantum Hadron Dynamics (QHD) model @, 0]. The 8 meson effects are also 
investigated in other models, such as a chiral SU(3) model Q, a relativistic point coupling model Q, relativistic 
transport model [l(| and so on. But there is no similar work in the quark-meson coupling (QMC) model yet, so we 
will investigate the 8 meson effects by using this model in this paper. 

The quark-meson coupling model was proposed by Guichon in 1988 [ll[ where nuclear matter is described as 
nonoverlapping MIT bags interacting through the exchange of mesons in the mean-field approximation. The effective 
nucleon masses in the QMC model are obtained self-consistently at the quark level, which is an important difference 
from QHD model. The model is refined by including nucleon Fermi motion and center of mass corrections to the 
bag energy by Fleck (l^ |. Jin and Jennings introduced the density-dependent bag constant, which is called modified 
quark- meson coupling (MQMC) model, to get larger scalar and vector potentials compatible with experiments [HI]. 
Furthermore, the MQMC model possibly includes the effects of quark-quark correlations associated with overlapping 
bags which was missing in the original QMC model, therefore it is applicable at the densities appropriate to neutron 
stars. The (cr*,4>) meson fields are incorporated to account for the strong attractive AA interaction observed in 
hypernuclei which cannot be reproduced by the (er, lj, p) only in MQMC model The MQMC model gives a 

satisfactory description of finite nuclei [l5| and nuclear matter [l6| . and it is widely used in nuclear physics. For 
example, the temperature effects of nuclear matter pjj , K condensation [lj| EH HH > trapped neutrinos [2(| , strong 
magnetic field [2l[ and deconfined phenomena [12] in neutron stars are all investigated in the MQMC model. 

In this paper, we extend MQMC model to incorporate 8 meson field, in which the density-dependent couplings 
between baryons and scalar mesons are calculated self-consistently. The model parameters are determined by the 
properties of symmetric nuclear matter and pure neutron matter. Then the influences of leptons, baryons and (tr*, 4>) 
mesons on the 8 meson effects are discussed. 



2. THE MODEL 

The modified quark-meson coupling model is extended to include the 8 meson field. 8 meson couples only to u 
and d quarks, because it is built out of nonstrange quarks, a* and <f> mesons are also incorporated which couple only 



2 



to the s quark in a hyperon bag. So there are isoscalar scalar mesons er and a*, isoscalar vector mesons u> and <f>, 
isovector scalar meson 6 and isovector vector meson p in our present model. 

In the mean field approximation the Dirac equation for a quark field of flavor q = (u, d, s) in the bag for the 
hadron species B = (p, n, A, £ + , £ , £~~, H , E - ) is then given by 



vy ■ d - (to, - g q a a - g%, a* - g q s h q h) ~ 7° (pfao + 9%<\>v + 9 q h q Po3j 4> q B (r, t) 



= 0. 



(1) 



Here I^ q is the isospin projection of quark q; g$ is the coupling constant between quark q and 5 meson, S3 denotes 
expectation value of the isospin 3rd-component of S meson field, and the other symbols are the same as in [23j |. The 
normalized ground state is solved as 
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with Rb is the bag radius of baryon B and x qB is the dimensionless quark momentum which can be determined by 
the linear boundary condition 



The effective quark mass is 



jo(x qB ) = f3 qB ji(x qB ) 
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(G) 
(7) 



The energy of a MIT bag for baryon B is then given by 

^ n qB Vt qB - z B 



pbag _ q 
B ~ 



Rb 



~-kR%B b (a,a*,S 3 ) 
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where n qB is the number of constituent quark q in baryon B, z B is the zero-point motion parameter and B B is the 
medium dependent bag parameter. The ansatz for the coupling of bag parameter to the scalar fields a, a* [241 ] is 
extended to 63 



B B (er, a* , 63) = B a exp \ - jj- 



bag * 1 \ ^ f bag 1 bag 7- c 
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where Bq is the bag constant in free space, M B is the bare mass of the baryon B, and g^ as , g^* 6 and gj? ag are real 
parameters. 

After the corrections of spurious center of mass motion, the effective baryon mass is given by 



where 



l qBX q B 



(10) 
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The bag radius Rb could be obtained through the minimization of the baryon mass with respect to the bag radius 



dM* E 
~dR~n 



= 



(12) 



Consider an many-particle system consisting of the full baryon octet which interact via a, a* , uj, <fr, 5, p meson fields. 
The Lagrangian density is 



ij^ - M* B (a, a*,5 3 ) - ( Q^yTf + 9 p b^- ■ P^ + 9^^ 



-d u 5-d^5 + d u a*d^a' 
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(13) 



where / = (e,p). Then from Eq. (jTUJ) and (H2J), we can derive the equations of the motion for the meson fields in 
uniform static matter: 
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Here k B is the Fermi momentum of the baryon species B. The factors Cb(ct), Cb(<t*), Cb{^) are: 
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The scalar density of quark q in the bag B are 

_ n qB /2 + R B m q (n q B-i) 
qB ~ n qB {n qB - 1) + R B m*/2' 



q = (u,d,s), 



At last, there are two conditions left: 



charge neutrality: ^ <lBk% = fc, 



13 equilibrium: p B = Mn - 9sMe, Mm = Me 



(20) 
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where qs and ps correspond to the electric charge and chemical potential of baryon B, respectively. The energy 
eigenvalue of Dirac equation for baryon B and lepton I are: 

e B = \Jk 2 B + + g^BUJo + g^B^o + g P BhBP03, (26) 

e, = ^/kf+mf (27) 

Then the Fermi momentum can be obtained from the equations 

e B (k B ) = p B (28) 

ei{h) = W (29) 

After the meson fields (c, a* , to, <f>, 5 3 , P03), Fermi momenta (k B , k{) and effective masses M B are obtained by solving 
the Eqs. (fl^ - fTS")) , (|28p - (l2"9")) and (jT3J) self-consistently at a given baryon number density 



= iE 6 ^ (30) 
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where b B is the baryon number of baryon B, we can obtain the total energy density and pressure: 
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3. PARAMETERS AND CALCULATION DETAILS 

Take the current quark mass to be m u = rrid = and m s — 150 MeV. Small current quark mass for the non-strange 
flavors or other values for the strange flavor lead only to small numerical refinements [191 ]. The meson masses are 
m a = 550 MeV, m* = 980 MeV, m p = 775 MeV, m s = 985 MeV, m u = 783 MeV, = 1020 MeV, respectively. 

Assume a, uj, p, S mesons couple only to the it, d quarks and cr*, <f> mesons couple only to the s quark, we have 

9<r = 9u = 9s = 9 P = 9a- = 9a* = 94, = 9<t, = (33) 
By assuming the SU(6) symmetry of the simple quark modelflij 

9 d a = <£, 9%* = V2 9 a, g d s = gl, gt = g%, g% = V2g%, g* = g u p (34) 

we can get the relations 

= ^g U A = ^g^s = g^s = gZ (35) 

g<t>A = g<f,x = -fl^s = V2gZ, g<t>N = (36) 



9 p n = g P A = g P T, = g P s = g u p (37) 



To reduce parameters we set 



bag bag bap 
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(38) 



The free nucleon zero-point motion parameter znq and the free bag constant Bq are fixed to reproduce the free mass 
of nucleon myv = 939 MeV with the minimization condition (fl"2"|) at a free bag radius Rno = 0.6fm. Then the free 
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TABLE I: The zero-point motion parameters zbo and bag radii Rbo in free space are obtained to reproduce the free space mass 
spectrum after the parameters — 188.102 MeV and zno = 2.030 have been fixed by the properties of nucleon. 





Mb (MeV) 


ZBO 


R B0 (fm) 


A 


1115.68 


1.815 


0.643 


£+ 


1189.37 


1.638 


0.669 


E° 


1192.64 


1.630 


0.670 


£~ 


1197.45 


1.612 


0.672 


CO 


1314.83 


1.501 


0.689 




1321.31 


1.483 


0.689 



TABLE II: Four independent coupling constants are fixed to reproduce the symmetric nuclear matter binding energy B/A — 16 
MeV, symmetry energy a sy m = 32.5 MeV and compressibility K = 289 MeV at saturation density po=0.17 fm~ 3 . 8 meson 
coupling constant is set and 4.2. 
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9 U P 







0.980 


2.705 


7.948 


2.278 


4.2 


0.980 


2.705 


10.217 


2.278 



zero-point motion parameters zbo and free radii Rbo of other baryons are obtained by reproducing the free baryon 
mass Mb with the minimization condition (|12|) . They are all listed in the Table [0 

Four independent coupling constants g^iff^'Sp and g^ ag can be adjusted by reproducing the symmetric nuclear 
matter binding energy B/A=1Q MeV, symmetry energy a sym =32.5 MeV and compressibility K = 289 MeV at 
saturation density po = 0.17 fm~ 3 , as listed in Table ILT1 The S meson-quark coupling constant is constrained in an 
range of < g$ < 4.2 so that the pure neutron matter EOS is consistent with the experimental flow data in heavy-ion 
collision [25|, which is shown in the upper panel of Figure. [TJ the EOS for symmetric nuclear matter is also shown in 
the lower panel and we can see that it is also consistent with the experimental flow data in heavy-ion collision. 

The equilibrium properties of neutron stars are obtained by solving Tolman-Oppenheimer-Volkoff (TOV) [26J equa- 
tions 

dP(r) = G[e(r) + P(r)] [M(r) + 4^r 3 F(r)] 
dr r 2 [l-2G¥(r)/r] 1 ' 

= We(r). (40) 

dr 

The Baym-Pethick-Sutherland model [27j is used to describe the EOS at subnuclear densities. 



4. RESULTS AND DISCUSSIONS 



Four cases in Table HTT1 are investigated: (1) pure neutron matter denoted by nn; (2) /3-equilibrium nucleonic matter 
denoted by np; (3) /3-equilibrium hadronic matter composed of baryon octet without (a*, <fi) meson fields, denoted by 
npH; (4) The same as in Case (3) with two additional meson fields (a*, 0), denoted by npH*. 

The meson fields for npH* are shown in the left panel of Figure. [2] We can see that 8 meson field decreases (a, u>) 
fields while increases (a*, (f>) fields. This is because the S meson increases the strange number in nuclear matter, which 
is shown in the right panel of Figure. [U and (a* , (f>) couple only to s quark. The 6 meson increases po3 meson field, 
and the effect becomes smaller when the S meson field decreases as baryon density increases. 



TABLE III: The cases we study in the paper. H represents hyperons (A, E + , E°, E , H°, 3 ). 
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a,u,p(,5) 


a,u,p(,S),a*,<j> 



G 
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FIG. 1: (Color online) The EOSes obtained in MQMC model for pure neutron matter and symmetric nuclear matter. The 
upper magenta hatched area and the lower orange hatched area correspond to the pressure regions for neutron matter after 
inclusion of the pressure from asymmetry term with strong density dependence and for symmetric nuclear matter consistent 
with the experimental flow data, respectively [25j]. 




FIG. 2: (Color online) The left panel are meson fields as functions of baryon density for npH* . The right panel is the ratio of 
s quark to total quark in nuclear matter versus baryon density. The upper two curves are for npH and the lower two are for 
npH* . 



Let's look at the compositions of nuclear matter for npH* in Figure. [3J 5 meson decreases the neutron fraction 
while increases the proton and lepton fractions when p > pq. From the right panel of Figure. Uwe see that S meson 
decreases the effective mass of neutron, which makes the neutron fraction fall when the density exceed some critical 
density which is approximately nuclear matter density po as shown in Figure. [3l The proton fraction goes up because 
the similar reason, and the charge neutrality condition requires larger lepton fractions. 

E _ ,S _ (negative isospin projection) and A, E° (zero isospin projection) appear earlier when S meson are included, 
but the appearance of E + (positive isospin projection) is postponed. From equations (j!?5j) . (f^H|) and we know 
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FIG. 3: The compositions of neutron stars in npH* as a function of baryon density and the vertical lines represent the central 
baryon densities of neutron stars with maximum mass. 
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FIG. 4: (Color online) The left panel is the changes of chemical potentials of neutron and lepton and the quark-meson exchange 
potentials versus the baryon density p when the 8 meson are included. The right panel is the differences of baryon effective 
masses. Both are for npH* . 



that the fraction for baryon B is determined by (u n — qBPe), (g^B^o + g^B^o + gpBhBPm) and M Bl which are all 
shown in the Figure. |H We see that the changes of g^o and g^<p are proximately offseted; M B and g^pm (Compare 
Figure. |4] with Figure. [21 we can see that the change of g p po3 mainly origins in the large change of quark-p meson 
coupling constant g™) change obviously. The changes of g P Bl3BP03 and M B are isospin-dependent, so whether the 
hyperon appears earlier is determined by its isospin projection. The critical density of its appearance only shifts a 
little except 3 - , since the changes of g P BhBPoz and M B are almost the same. 

The S _ population with S meson is larger than that without S meson at the beginning, but it becomes smaller 
than that without S meson because of the appearance of S - . The reasons are that charge neutrality can be kept more 
economically by the larger mass particles with the same charge, and the S meson decreases Mt more than M~_ (right 
panel of Figured]) since the isospin projection of S~ is —1/2 and £~ is —1. S meson increases S~ population obviously 
larger than other hyperons since it decreases £~ population. The appearance of S _ makes the lepton fraction begins 
to fall, which can also be explained by charge neutrality condition. There is another interesting phenomenon that 
the E + may not appear in neutron stars with S meson while its fraction could exceed 1% for neutron stars at the 
maximum masses without S meson. 

The EOSes for nn, np, npH, npH* are plotted in Figure. [5] The effects of 6 meson can be seen clearly from this 
figure: The S meson makes the EOS of nn stiffer similar as in QHD model (281 ) - For np, the S meson stiffens the EOS 
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FIG. 5: (Color online) The equation of state, pressure P versus energy density e for all cases we study in this paper. The 
causual limit (P — e) is also shown. 
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FIG. 6: (color online) The mass-radius relation for neutron stars in nn, np, npH*, npH. The filled triangles, inverted triangles, 
squares and circles represent the positions of maximum masses for nn, np, npH* , npH EOSes, respectively. A gravitational 
redshift of 0.35 [l^] (observed from EXO 0748-676) is shown in this figure. Mass-radius constraints from thermal radiation of 
isolated NS RX J1856. 5-3754 ^ (cyan hatched area) and from QPOs in the LMXBs 4U 0614+09 [H| (magenta hatched area) 
and 4U 1636-536 [32| (orange hatched area) shall be regarded as separate constraints to the EOSes. The mass of PSR 1913+16 
[33| is also shown in this figure. 



at low density while softens at high density. The density-dependent coupling constants are introduced additionally in 
QHD model [7[ to get the similar results, but the density-dependence of couplings between scalar mesons and baryons 
are obtained self-consistently in our paper. If hyperons are taken into account, the EOSes with S meson suffer a 
transition to nucleon-hyperon phase at some density and become softer, this can be seen from the EOSes of npH and 
npH* in Figure. [5] clearly. (ct*,</>) mesons obviously stiffen the EOSes as in Ref . (lij . but their influences on the S 
meson effect could be neglected. Since (a*,<j>) meson fields couple only to s quark and S meson couple only to (u, d) 
quarks, 5 meson has no direct influence on (a* , <fi) meson fields, which can also be seen from the lower panel in Figure. 
[5] as mentioned above. 

We find that no matter whether hyperons are positive, negative or neutral, their inclusions can make the EOSes 
with 5 meson become softer. This result probably reveals that it is the strange quark makes the EOSes with 8 meson 
become softer. That is to say that it is the strange quarks in hyperons results in reversed direction changes of EOSes 
if compared with a nucleonic star. 

The mass-radius relation of neutron stars are shown in Figure. [5] We see that the 5 meson increases the maximum 
masses of neutron stars for all cases we studied in this paper. This is different from QHD model in which the maximum 
mass decreases for np with density-dependent couplings Q and npH when 5 meson are included 28]. The S meson 
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TABLE IV: The maximum masses of neutron stars and the corresponding radii -RAf max , central baryon density p c , central 
energy density e c , central pressure P c for different EOSes. 





gs 


M max /M s i?M max 
(km) 


Pc/P0 


(MeV for 


Pc 

- 3 ) (MeV fin 


nn 


4.2 


2.275 


12.07 


5.1 


1093 


418 


nn 





2.147 


11.50 


5.6 


1213 


459 


np 


4.2 


2.045 


11.30 


6.0 


1274 


446 


np 





2.012 


10.95 


6.3 


1352 


493 


npH* 


4.2 


1.556 


11.56 


6.0 


1157 


201 


npH* 





1.543 


11.10 


6.5 


1274 


236 


npH 


4.2 


1.509 


12.03 


5.2 


986 


140 


npH 





1.491 


11.54 


5.8 


1106 


166 



enlarge the radii of neutron stars about 0.5 km for stars with M > M s , this is an obvious change considering the same 
EOS at low density are used for all cases. Another conclusion is that the central density of neutron star becomes 
about 0.5 po smaller when 6 meson is included. These can be seen from Table HVl clearly. Some observation values are 
also displayed in Figure. [6] We can see that all cases are compatible with the observations from PSR 1913+16 [33| 
and 4U 0614+09 [3l[, but npH* and npH neutron stars might be ruled out by neutron star 4U 1636-536 [HJ or EXO 
0748-676 [Hj]. To show quantitatively the 5 meson effects on neutron stars properties, the maximum mass M max and 
the corresponding radii i?M max , central baryon density p c , central energy density e c , central pressure P c are listed in 
Table ITVl for all cases. 

5. CONCLUSIONS 

We have investigated the 6 meson effects on neutron star proporties within the modified quark-meson coupling 
model. We sum up the conclusions in four aspects: 

(1) 6 ^strangeness: the 6 meson can make the pure nucleonic matter EOSes stiffer, while make the hyperon matter 
EOSes softer and this could be explained by the appearance of strange quarks in hyperons. 

(2) 6 ^-deptons: the 6 meson results in opposite effects on the EOSes of /3-equilibrium nuclear matter such as np, 
npH and npH* compared with the EOS of nn, which is similar to the influence of hyperons, but the effect is smaller. 
This is because of density-dependence of the couplings between baryons and scalar mesons. 

(3) 6 «-» (a-*,(j>): 6 meson has no direct influence on s quark, so it has little effect on (<j*,(p) meson fields. As a 
result, (<j*,<t>) mesons have no obvious influence to the 6 meson effect on EOSes, although (<J*,<j>) stiffen the EOSes 
obviously. 

(4) <5 <-> neutron star properties: the 6 meson can increase the maximum masses of neutron stars, decrease the 
corresponding both baryon density and energy density. The radii become about 0.5 km larger for stars with M > M s . 
It changes compositions of neutron stars: decrease the neutron fraction and increase the proton and leptons fractions 
when p > po! make the abundance of S _ ,S° and S + larger, while the abundance of S _ smaller; and increase the 
strange number of neutron stars. It can also make the isospin dependent physical quantities splitting, such as effective 
baryons masses. 
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